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INTRODUCTION 


Cryogenic wind tunnels have become important tools for high Reynolds number 
research. (See ref. 1 by Kilgore et al. ) A cryogenic wind tunnel (T'2) has been 
built at the Toulouse Research Center, Toulouse, France. (See ref. 2 by Blanchard, 
Dor, and Breil. ) The 0.3-Meter Transonic Cryogenic Tunnel (0.3-m TCT) , built as a 
pilot facility, is at the langley Research Center. The National Transonic Facility 
(NTF) at the Langley Research Center (see ref. 3 by Fuller) , which is ejipected to 
become operational in 1983, is a closed-circuit, fan-driven, continuous-flow, pres- 
surized cryogenic tunnel operating at pressures up to 9 atm (1 atm = 101.1 kPa) and 
having a 6.25-m^-area slotted test section that is 7.6 m in length. The facility 
will operate at tenperatures ranging from ambient down to about 100 K. Liquid nitro- 
gen (LN 2 ) sprayed into the tunnel upstream from the fan nacelle performs three 
cryogenic-cooling functions: (1) initial cooldown, (2) stea<ty-flow tenperature regu- 

lation by balancing the energy added by the fan, and (3) primary tenperature control 
in passing from one steady-flow condition to another. Total-pressure control is 
accomplished by venting tunnel gaseous nitrogen (GN 2 ) to the atmosphere. Both LN 2 
and GN 2 flows are regulated by means of servo-control valves. Fan-motor speed deter- 
mines coarse Mach number variation, whereas fan inlet guide vanes furnish fine ver- 
nier Mach number variation. The tunnel shell has interior thermal insulation with an 
aluminum inner liner. This limits atmospheric-heat penetration to a negligible 
amount and reduces energy consumption. See figure 1 for a schematic diagram of the 
NTF. 


Inlet guide 



The high cost of liquid nitrogen makes efficient operation of a cryogenic tunnel 
inportant. In order to minimize operating cost, it is necessary that such a facility 
be equipped with automatic controls which perform the following functions: 

(1) Maintain stea<^ flow at the most efficient setting during test dwell. 



(2) Drive the system along the most efficient transition path from one stea<fy- 
flow setting to the next as quickly as possible. 

Modern control-system synthesis techniques require that the dynamic process to 
be controlled be modeled analytically either in a discrete formulation as a system of 
difference equations or in a continuous formulation as a system of differential equa- 
tions in state variable form. Hie thermodynamic and fluid-dynamic processes in the 
wind tunnel are described by the three-dimensional Navier-Stokes equations which 
require numerical solution. Hie complexity of this three-dimensional formulation is 
beyond the scope of modern distributed-parameter control synthesis techniques, and 
numerical solution of the equations requires lengthy, expensive computation. Conse- 
quently, the control- system designer must resort to a simplified approach, such as 
segmentation of the process into a set of lunped parameter stations which may then be 
treated as a high-order multivariable system with time delays. Because the accuracy 
of such an approximation is limited, a need exists for a dynamic model whose accuracy 
is adequate for control- system design and performance evaluation, but whose numerical 
solution will execute on a digital conputer at reasonable speed and cost. The one- 
dimensional distributed-parameter dynamic model described herein is intended to 
satisfy this need. 


SYMBOLS 


^S 

Ap 

a 

B 

b 

Cd 

“^DMX 

‘^ORF 


'XM 


'M 


diffuser area ratio 

2 

pro jected-strut frontal area, m 

test-section cross-sectional area at throat, m^ 

local speed of sound, ra/sec 

body force vector in three-dimensional space, N 

_ 1 


2y 


(see eq. (A1 4) ) 


drag coefficient 

variable model-blockage flow coefficient 
flow coefficient for reentry mass flow 

slot (plenimi test-section) pressure-difference flow coefficient, 
kg-m^/N-sec 

heat-transfer correction factor 

variable defined in equation (93) 

artificial- viscosity coefficient 

boundary-layer flow coefficient 

specific heat of metal liner, kJ/kg-K 

specific heat of gas at constant pressure, kJ/kg-K 
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V,n 


u ,n 


E 


PL 


f 

G 

GN2 

H 

h 

j' 

K 

^6 

k 

’"h 

L 

LN2 

M 


m 

m' 

N 


Npr 

^Re 


specific heat of gas at constant volume, kJ/kg-K 
vector difference defined in equation (95) 
tube diameter, m 

difference defined in equation (94) 
plenum internal energy, kJ 
internal energy, kJ/kg 
vector defined in equation (76) 

surface force vector in three-dimensional space, N 
function 

vector-function relation of fan 

gaseous nitrogen 

vector defined in equation (77) 

enthalpy, kJ/kg 

total enthalpy of LN 2 , kJ/kg 

momentum-flow rate per unit length, kg/sec 

diffuser loss factor 

diffuser correction factor defined in equation (37) 

constant defined in equation (45) 

thermal conductivity of gas, kJ-m/K-sec 

diffuser correction factor defined in equation (34) 

length, m 

liquid nitrogen 

Mach number 

mass, kg 

mass-flow rate per unit length, kg/m-sec 
number of quantized stations 
Prandtl number 
Reynolds number 



V,n 


P 

K 

Q 

q 

q' 

R 

% 

Rh 


M 


T 

T. 

t 

U 


U, 


T 


U 


V 


V 


vector defined in equation (96) 
static pressure, atm (1 atm = 101.1 kPa) 
corrected total pressure 

external heat flow into control volume, kJ 
heat-flow rate, kJ/sec 

heat-flow rate per unit length, kJ/m-sec 

heat transfer from liner, kJ/m-sec 

gas constant, kJ/kg-K 

constant defined in equation (23) 

wall resistance, km /kJ 

diffuser-inlet equivalent radius, m 

fan conpression ratio 

diffuser pressure ratio defined in equation (33) 
corrected diffuser pressure ratio 

2 

surface vector in three-dimensional space, m 

entropy, kJ/K 

gas temperature , K 

liner temperature, K 

time, sec 

liner thickness, m 
velocity vector 

3 

heat-transfer coefficient, kJ/km 
artificial-viscosity vector 
velocity, m/sec 

state vector defined in equation (73) 

value of state vector V at station x^ and time tj^ 
specific volume, m /kg 
volume, m^ 
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w 

w 

''dm 

''f 

''G 

''n 

''ss 

WrpH 

X 

a 

Y 
A 
6 

6 * 

e 

\ 

29 

V 

P 

Pm 

a 

(jO 


external work done by control volume, kJ 
mass-flow rate, kg /sec 

boundary-layer conponent of slot-flow rate, kg/sec 

model -blockage conponent of slot -flow rate, kg/sec 

extracted conponent of reentry mass-flow rate, kg/sec 

GN 2 ^®rit flow rate, kg/sec 

LN 2 -injection flow rate, kg/sec 

supersonic conponent of slot-flow rate, kg/sec 

average test-section mass-flow rate, kg/sec 

spatial coordinate, m 

inlet-guide-vane position, deg 

ratio of specific heats 

differencing increment 

increment used in equations (33) and (39) 
boundary-layer displacement thickness, m 
increment used in equations (38) and (39) 
fan efficiency 

total diffuser expansion angle, rad 
boundary-layer momentum thickness, m 

p 

coefficient of viscosity, N-sec/m 
fan speed, rpm 
density, kg/m^ 
liner density, kg/m^ 

p 

normal stress, N/m 

p 

cross-sectional area, m 


Subscripts: 


DM model blockage 


E slot exit (plenum test section) 
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EX test-section exit 

E fan outlet 

G GN 2 (gaseous nitrogen) 

N LN 2 (liquid nitrogen); tunnel outlet, for example, in equation (74) 

ORF plenum test section 

PL plenum 

R flow loss 

RE reentry 

s smoothing 

TS test section 

t total value 

XX x-component of stress of force 

Superscripts; 

- average value 

~ predicted value 

■' corrected value 

Special symbol; 

V vector del operator 

A dot over a symbol denotes a derivative with respect to time. 


DEVELOPMENT OF ONE-DIMENSIONAL MODEL 
General Assumptions 

The wind tunnel is modeled as a one-dimensional tube of varying cross-sectional 
area, in which the flow is assumed to be uniform across every cross section. Con- 
sequently, rotational effects and mixing caused by turns in the tunnel are neglected. 
In addition, longitudinal mixing caused by diffusion and turbulence is not modeled. 
Although viscous shearing stresses are neglected, frictional momentum losses at the 
walls are included. 

Thermodynamic properties are computed by means of ideal-gas laws. Real-gas 
effects, important at low temperatures, could be readily included. Heat transfer is 
assximed to occur only between the gas and the tunnel inner liner; heat penetration 
from the outer shell through the insulation to the liner is neglected. Liquid- 
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nitrogen evaporative effects on total pressure, which are observed primarily at the 
lowest temperatures, have not been considered. 

The tunnel fan is represented as an algebraic functional relationship with no 
inherent flow dynamics. Thus, the fan model is instantaneous, exhibiting no time 
delay. Inlet-guide-vane actuator dynamics are included. The fan model determines 
the functional relationships between tunnel inlet and outlet flows, which effectively 
closes the tunnel circuit. 

The plenum is modeled as a lumped voliime attached to the distributed test 
section. All external processes, including actuators and automatic controllers, are 
modeled as lumped-parameter dynamical systems. The nimierical solution of the set of 
ordinary differential equations which describe the lunped-parameter processes is 
synchronized in time with the numerical solution of the partial differential equa- 
tions of flow. This synchronized combination of one- and two-dimensional systems is 
well-suited for the CDC^ CYBER 203 vector-processing digital computer (ref. 4) on 
which the solutions are obtained. 


Derivation of One-Dimensional Flow Equations 

The three-dimensional equations of fluid flow in integral form are applied to a 
differential element of a one-dimensional tube of varying cross-sectional area in 
order to derive a one-dimensional system of partial differential equations analogous 
to the three-dimensional Navier-Stokes equations. The three basic conservation laws, 
expressed in vector notation (ref. 5), are given as follows: 

( 1 ) Continuity : 

J pU • dS + J p du = 0 (1) 


( 2 ) Linear momentum : 


I 


pU 


U • ds + 


9t 



dy 



B du 


(2) 


(3) Energy: 



dS + 


5t 



di; 


dQ dW 
dt dt 


where 

p density of gas 

U velocity vector of mass flow 


(3) 


^CDC: 


Registered trademark of Control Data Corporation. 
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S surface vector of volvune 

V volume 

Fg surface force vector 

B body force vector 

p static pressure 

e^ total internal energy of gas 

Q external heat flow into control volume 

W external work done by control volume 

By referring to figure 2, consider an incremental element of tube of length Ax 
with flow entering at x^ and exiting at X 2 * 



Figure 2.- Elemental flow volume u. 


The inviscid one-dimensional form of equations (1), (2), and (3) (ref. 6) is, 
respectively , 



+ 


9x 


( puu) 


= 0 


(4) 


dt 


{ puco) 


^ / 2 
. ^(pu » 


+ pw) 


= P 


dco 

dx 


(5) 
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and 




( 6 ) 


where w is the cross-sectional area of the tube and u is stream velocity. 

Equations (4), (5)/ and (6) are augmented to account for external mass, 

momentum, and heat transfer, respectively, into the element of volxime. Let m' , 
j', and q' denote mass-flow rate, momentum-flow rate, and heat-flow rate per unit 
length, respectively, into the element of volvmie. Hie surface integral in equa- 
tion (1) for conservation of mass with mass flow into the element of volume becomes 

^ pu • dS = Ax - m' Ax (7) 


so that equation (4) becomes 




+ ^(puw) = m' 


( 8 ) 


Similarly, the addition of momentum flow j' Ax and heat flow q' Ax into the 
element of volume results, respectively, in the equation for conservation of momentum 


-g^(puo)) 


5 , 2 

^(pu 0. 


+ pco) = P 


do) 

dx 


(9) 


and in the equation for conservation of energy 


at 




pco) + ^(e^puto + pu(o) 


q* 


( 10 ) 


Equations (8), (9), and (10) are in the form employed by Carriers (ref. 7) in a 

dynamic analysis of an injector-driven wind tunnel that was solved by means of the 
method of characteristics. 

In the present study, it is desired to account for heat transfer within the gas 
due to gas thermal conductivity. From Fourier's law of heat conduction, it follows 
that the rate of heat transfer through surface S into volume u is 

I k VT • dS (11) 

where k is the thermal conductivity and denotes the component of heat transfer 

due to gas thermal conductivity. Because gas-to-gas heat transfer does not occur 
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along the outer surface of the surface integral, equation (11) is evaluated only at 
surfaces and CO 2 (fig* 2) so that 


dt 


k 


k(02 


dT 

5x 



( 12 ) 


(10)) then becomes 

= q' (13) 

In this analysis, it is desired to consider the effects of fluid viscosity. The 
normal stress for one-dimensional flow is (ref. 8) 


The equation for conservation of energy (eq. 




putjj + puo) - ko) 


^(e^pw) 


a 

XX 




(14) 


where p is the coefficient of viscosity. In the present one-dimensional approxi- 
mation, shearing stresses are neglected. The surface force at points 

and X 2 (see fig. 2) due to normal viscous forces is of the form 


F 

XX 


too 


XX 



(15) 


Consequently, the equation for conservation of momentum (eq. (9)) then becomes 

■^(puco) t + pu) - i pco = P + j. ( 

Empirical momentum losses obtained from data reported by Rao (ref. 9) will be 
introduced through the j ' term. 

The rate of flow work done at surfaces and u)^ by the normal force F 

. . ^ . 1 2 XX 

in equation (15) is 


dW 


dt 


w 


= F • U| 

XX 


4 9u 

w = 3 ^ 


(17) 


The energy equation for viscous flow with heat addition 


and thermal conduction is 




e^puo) 


+ puo) - ko) - 4 
dx 3 


puw 


5u 
9x , 


9t 


(e^pw) 


q’ 


(18) 
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Equations (8), (16), and (18) cortprise the conplete one-dimensional partial dif- 
ferential equations of flow that approximate the three-dimensional Navier-Stokes 
equations . 


Heat Transfer Between Gas and Tunnel Wall 

The one -dimensional time -dependent equations of flow developed in the preceding 
section include gas-to-gas heat-transfer effects. It will be found that heat trans- 
fer between the gas and the shell of the tube, corresponding to the wind-tunnel metal 
liner, is of primary importance. Let the rate of heat transfer from the liner to the 
gas (in joules per second per unit length) be denoted by qg, the metal tenperature 
by Tj^, the gas tenperature by T, and the heat-transfer coefficient by U,p. For an 
annulus of tube of diameter d, the heat-transfer rate is 


q' = -n:d(T^ - T)U^ 
S M T 


(19) 


An empirical relationship for U,j given by McAdams (ref. 10) is 


1 



1 


„ „ k 0.8 
0. 026 - 

d Re 


N 


.0.4 

Pr 




( 20 ) 


where 

k thermal conductivity of gas 

Njjg Reynolds number 

Np^ Prandtl number 

R^ wall resistance 

It has been found that a heat-transfer correction factor Cp is necessary to account 
for the additional surface area in the nacelle and turning-vane sections. The factor 
is employed as a multiplier so that equation (19) becomes 


q^ = TCd(TM - T)U^Cj^ (21) 

Consider next the tunnel inner-liner temperature dynamics. Heat transfer from 
the liner to the gas occurs in accordance with equation (21). lateral heat transfer 
within the liner from warmer to cooler regions is neglected. Consequently, the rate 
of change of the liner tenperature T is 


T = 
M 


(T - T^)U^ 


( 22 ) 


1 1 



where 


Cjj is the specific heat of the liner material and 


where Pj^ is the liner density and tj^ is the liner thickness. Equation (22) is 
employed as the governing equation for Tjj. 


ESTIMATION OF FLOW LOSSES 
Introduction of Losses Into Equations of Flow 

In the one-dimensional equations of flow, the viscous shear effects at the 
tunnel walls are neglected. However, frictional effects, especially in the test 
section and diffuser section of the wind tunnel, are significant and must be 
accounted for. According to Rao (ref. 9), nearly 60 percent of the total aerodynamic 
loss at Mach 1 occurs in these two sections. Since, of course, energy is conserved 
in the test section and diffuser, the "total energy loss" to which Rao refers is a 
loss in available energy, which is manifested as an increase in entropy. The entropy 
gradient at steady flow ds/dx, obtained in the appendix in equation (A23), is given 
by 



Because of the assiamption that the tunnel liner is perfectly insulated, external 
heat transfer q* equals zero for steady flow. Furthermore, mass transfer term m' 
is zero everywhere except for GN 2 and LN 2 transfer. Therefore, the increase in 
entropy must be due entirely to a momentxim loss j' so that 


j ' < 0 


(25) 


The steady-flow expression for the gradient of total pressure given in equation (A18) 
is 





um' + (y - 1 )m‘ 


(j' 


(26) 


With m* 


and q’ both zero, equation (26) becomes 




dx 






( 27 ) 
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and because j' is negative, dp^/dx is negative. Rao (ref. 9) furnishes empirical 
relations for determination of total-pressure losses in the test section, in diffuser 
sections, at turning vanes, at cooling coils, and in the strut section of the NTF 
wind tunneli With Rao's relations, dp^/dx can be determined, and j' may be 
obtained from equation (27) for each of the tunnel sections as 



(28) 


where (dp^/dx)j^ is the total-pressure gradient due to flow losses obtained from 
Rao's relations. This value of j' is entered into the dynamic equations for 
unsteady flow, and dynamic effects in j* are neglected. 


Calculation of Diffuser Losses 

Total-pressure loss Ap^ along the length of a diffuser is given enpirically by 
the equation (ref. 9) 


- ft.l - - Pi 


(29) 


where 


K 

Pt,1 

Pi 

Pt,2 


diffuser area ratio 
diffuser loss factor 
total pressure at inlet 
static pressure at inlet 
total pressure at outlet 


Studies by Rao (ref. 9) and Henry et al. (ref. 11) show that for values of 
total diffuser expansion angle 29 lying between 4® and 10°, the diffuser loss 
factor K depends on the ratio of the inlet boundary-layer displacement thickness 
6* and the diffuser-inlet equivalent radius R^. For values of 6*/R^ less than 
0.030, K is well approximated by the relationship 


K = 0.075 + 5.88 


6 * 


(30) 


Substitution of equation (30) into equation (29) gives 


^t- <Pt,i -Pi>(' 


075 + 5.88 




(31) 
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Figure 3 shows a comparison of the diffuser pressure ratio 


= P 


t, 1^^t,2 


( 32 ) 


o 

cu 

i— 

cu 


Q 



Figure 3.- Comparison of computed and measured diffuser pressure loss in 
the Langley 0.3-Meter Transonic Cryogenic Tunnel. 


computed by using equation (31), with experimental values obtained by Rao (ref. 9) 
in the Langley 0.3-Meter Transonic Cryogenic Tunnel over the Mach number range from 
0.2 to 1.0. Agreement is good up to a Mach number of 0.8, above which the computed 
values of rj^ are too small. The ratio r^^, being nearly equal to 1.0, may be 
expressed by 


tj, • 1 + 6 


(33) 


where 6 is a 
data. Because 
tion factor to 


small increment and 
the values of r^^ 
be applied to rjj 


0 < 6 < 0.061, as computed from the experimental 
computed by equation (32) are too small, a correc- 
can be defined as 


k 


H 


Pt,2/P^,2 


(34) 
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where 2 corrected value of 2 ' giving a corrected diffuser pres- 

sure ratio 


r ' = k r = p /p ' 

H H H ^t,r^t,2 


(35) 


The corrected pressure difference can be expressed as 



P 


t f 


1 


(36) 


The ratio of Ap^ to Ap^ defines a correction factor Kjj. From equations (29), 
(32), (35), and (36), it follows that 


‘p; p;,2 - Pt.i _ - ’ . , ^ ^ 


' *pt ■ Pt,2 • Pt,i "''■'h’ - ' 


'h - ' 


Because is slightly greater than 1, let 




Substituting equations (33) and (38) into (37) gives 


(37) 


(38) 


for values of M > 0-6. Correction factor Kj^ equals 1.0 for M < 0.6. The com- 
putations of e, 6, and Kjj from the data of figure 3 are summarized in table I. A 
second— degree polynomial, fit to the values of Kjj obtained as a function of M 
for M > 0.6, is 


K « 1.0 + 1.375(M - 0.6)^ 
H 


(40) 


TABLE I.- CORRECTION FACTOR OF DIFFUSER PRESSURE RATIO 


M 


6 


e 

Kjj « 1.0 + ( e/6) 

Kjj » 1.0 + 1.375(M - 0.6)^ 

0.6 

1.025 

0.025 

1.0 

0 

1.00 

1.00 

.8 

1.042 

.042 

1.0029 

.0029 

1.07 

1.06 

.9 

1.052 

.052 

1.0076 

.0076 

1.15 

1.12 

1.0 

1.061 

.061 

1.0132 

.0132 

1.22 

1.22 
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which is tabulated in the final column of table I. Thus, is computed empiri- 

cally by the relation obtained from equations (31), (36), and (37) as 



where Kj^ is obtained from equation (40). 

The value of 6*/R^ is required as a function of Mach number M and Reynolds 
niamber Nj^^. Shapiro (ref. 12) obtains 6*/0^ as a function of M where 0^ is 
defined as the momentum thickness of the boundary layer entering the diffuser. A 
curve of 6*/0^ as a function of M, denoted by f 2 (M), is shown in figure 4 for 
values of M lying between 0 and 1.2. Rao (ref. 9) estimates 0^ by the relation 



Mach number 

Figure 4.- Variation of fj with Mach number. 



where L,j,g is the test-section length and is the drag coefficient for turbulent 

flow past an insulated flat plate. Shapiro (ref. 12) gives Cj^ as a function of 
Mach number and Reynolds number as 


C = 0.472(log.. + 0 . 2M^ (43) 

D 1 0 R© 

An ejqpression for 6* as a function of Mach number and Reynolds number can be 
obtained as 


6* = e, f^(M) = K- f_(M) 
12 o 3 


(44) 


where 


^6 


0.472(log^P 


^Re^ 


-2.58 



(45) 


and 


f^(M) = (1 + 0.2M^)"°‘'^^^ f2(M) (46) 

Values for £ 2 ^'^^ f 3 (M) are given in table II and appear as curves in fig- 

ures 4 and 5. A second-degree polynomial fit to f-j(M), denoted by f 4 (M) and 
tabulated in the rightmost column of table II, is 


f (M) = 1.286 + 0.305M 
4 


(47) 


The curves for f 3 (M) and f^(M) coincide in figure 5. 


TABLE II.- FLOW-LOSS RELATIONS AS FUNCTIONS OF MACH NUMBER 


M 

f2(M) 

f3(M) 

f4(M) 

0 

1.286 

U286 

1.286 

.2 

1.304 

U299 

U298 

.4 

1.358 

1.338 

1s335 

.6 

1.447 

1.401 

U396 

.8 

1.573 

1.487 

1.481 

1.0 

1.734 

U593 

1.591 

1.2 

U930 

1.715 

1.725 
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f = ^ (1.0 + 


CO 

CO 

o ( 


o / 

= 1.286 + 0.305M^ — 

y 

y 

y 


Mach number 


Figure 5.- Variation of and with Mach number. 


For the second- and third-leg diffusers and the rapid diffuser in the NTF, Rao 
(ref. 9) obtains a constant diffuser loss factor K of 0.32. In these locations, 
equation (31) simplifies to 


- 0.32P, (l . a^-L 


The portion of the fan nacelle beyond the fan location is an annular diffuser whose 
total-pressure loss is computed by Rao for a constant drag coefficient to be 




where Cp = 0.158. 

Losses at turning vanes, cooling coils, and screens are also obtained by using 
equation (49), where values of are given as follows: 


Location of loss 


Turning vanes 

0.027 

Cooling coils 

12.0 

Screens 

3.0 


Pressure loss at the strut, also obtained from equation (49), is weighted by the 
ratio of projected strut frontal area Ag to test-section cross-sectional area A,j,, 
thus giving 


Apt = 



(50) 


where = 0.027. Pressure losses due to the model and sting are confuted similarly 

as a function of angle of attack. 

In each aforementioned case, the total-pressure loss at each diffuser or con- 
stricting element is uniformly distributed along its length to obtain a value for 
(<^Pt/dx)j^ to be used in equation (28). 

SLOTTED-TEST-SECTION DYNAMICS 

The flow processes occurring within the slotted test section and plenum of a 
transonic wind tunnel are not understood well enough to permit construction of an 
accurate dynamic model. It is clear that the process is asymmetric and three- 
dimensional with sharp gradients and finely detailed flow patterns. An adequate 
time-dependent model would require a three-dimensional analysis with a finely spaced 
grid, which would require lengthy execution on a high-performance digital conputer 
for solution. Such an analysis has not been attenpted in the literature as of this 
writing. 

In this study, a one-dimensional approximation to slotted-test-section dynamics 
is adapted from a lumped model of the NTF slotted-test-section dynamics developed by 
Gumas. (See ref. 13.) The test section is divided into an exit region (from which 
slot flow exits into the plenum) which is followed by a reentry region (where flow 
reenters the test section from the plenum) . Slot exit flow is considered by Gumas 
to consist of three components caused by boundary layer, model blockage, and super- 
sonic effects. Let the lumped exit and reentry mass flows of Gumas be distributed 
along their corresponding regions. The plenum is considered to be a lumped volume 
with static plenum pressure uniformly distributed along the exit and reentry regions 
of the test section, in contrast to the treatment by Gumas of the test-section — 
plenum combination as a single lumped volume at uniform pressure. This separation of 
a luitped plenum from a distributed test section introduces another component of slot 
mass flow proportional to the difference in static pressure between test section and 
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plenum not included in the formulation of reference 13. The performance of the model 
for the subsonic case using these assunptions is satisfactory; however, above 
Mach 1.0 it is less successful, having a tendency toward numerical instability. 



Throat 


Strut 


Diffuser 


Pts 



A diagram of the throat, test section, and plenum regions appears in figure 6. 
Let Wg denote slot exit (test section to plenum) mass flow which, as mentioned 
before, according to the Gumas model consists of three components; 

(1) The boundary-layer conponent proportional to throat flow: 


w = C w 
B XM TH 


(51) 


(2) The model-blockage conponent, also proportional to throat flow and dependent on 
angle of attack; 


w = C w 
DM DMX TH 


(3) The supersonic component, when M,j,g > 1.0: 



(52) 


( 53 ) 
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where 


"TH 

Pt 

A.J, 

^XM 

^DMX 


average test-section mass-flow rate in throat region 
average test-section Mach number in throat region 
total pressure 
total temperature 

test-section cross-sectional area at throat 
boundary-layer flow coefficient 
variable model-blockage flow coefficient 


The additional pressure-difference flow component is computed as 


'^ORF “ ‘-ORF^PtS PpL^ 


( 54 ) 


where 


Pts 

PpL 

*^ORF 

Total 


average static pressure over exit-reentry portion of test section 
plenum static pressure 

slot (plenum test-section) pressure-difference flow coefficient 
slot-exit flow Wg is then 


’’’ '^SS '^ORF 


( 55 ) 


Note that if the plenum pressure is sufficiently greater than the test-section static 

pressure, w_ will be negative. Slot exit momentum flow j„ and heat-energy flow 

q_ are obtaxned from w_ as 
■‘‘E E 


fW U 

E TS 


^E 


(w > 0) 
E 

(w < 0) 

E 


( 56 ) 


where u^g is the average test-section stream velocity in the exit region and 




. h 
E TS 


Iw h 
E PL 


(w > 0) 

E 

(We < 0 ) 


( 57 ) 
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where 


h = G T 
TS p TS 


= c T__ 
PL p PL 


(58) 


(59) 


and Tpj^ is the plenum temperature. 

Reentry mass flow w^^, is computed from two components: (1) An extracted 

component proportional to test-section exit flow 


"'f = -""f^EX 


(60) 


where w^j^ is the average test-section mass flow in the reentry region; and (2) the 
pressure-difference flow component '^ORF* Hence, total reentry flow w^^ is 


RE 


- w 


ORF 


(61) 


Reentry momentum flow and energy flow are treated analogously in regard to slot-exit 
flow. The dynamic state equations for plenum mass and internal energy, respectively, 
are 


and 


where 


”*PL '^E ^ '^RE 


E = a + a 
PL ^E ^RE 


^RE^PL 




^RE 


'^reNs 


(w > 0) 
RE 


Plenum temperature is given by 


(62) 


(63) 


(64) 


T = E /c m 
PL PL'^ V PL 


( 65 ) 
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BOUNDARY CONDITIONS 


GN 2 Venting and LN 2 Injection 

Next, consider the external mass, momentum, and energy transfer caused by gas 
venting and liquid-nitrogen injection. Let the vented-gas mass-flow rate be denoted 
by Wq, and let the influent mass-flow rate per unit length be m^i. Assume uniform 
venting over length Lq so that m^ is given by 

m- = -w^/L^ (66) 


The influent momentum-flow rate per unit length jg is then 


jZ, = mZ.'i = “W^u/L„ 
■^G G G G 


(67) 


where u is the local stream velocity at the vent. 

Finally, the influent heat-flow rate per unit length is 


q' 


= m'h„ = -w^. 
G G G 




( 68 ) 


where hg, the local total enthalpy, is given by 


h 


G 


m . ■> 2 

c T + — u 

p 2 


(69) 


Corresponding equations for liquid-nitrogen injection can be obtained. Let w^ 
denote LN 2 -in jection flow rate (in kilograms per second) occurring uniformly over 
length L^^. Influent mass-flow rate per unit length denoted by m^ is then 

% = V-N 


Because the injected liquid nitrogen is assumed to have zero stream velocity at the 
point of injection, influent momentum-flow rate per unit length j^ is 




(71) 


Finally, influent energy per unit length is given by 




‘^N = %^N = -i;; 


(72) 


where h^ is the total enthalpy of liquid nitrogen. The value used for hj^ is 
-115.9 kJ/kg. Instantaneous vaporization and perfect mixing are assumed. 
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NTF-Fan Simulation 


A one -dimensional approximation to the NTF fan-flow process has been adapted by 
Gumas (ref. 13) from a three-dimensional fan model. In the approximation by Gumas, 
fan volumes are assumed to be insignificant relative to the tunnel volume; mass, 
momentum, and energy relations between fan-inlet flow and fan-outlet flow are assumed 
to be algebraic functions dependent on fan speed v and inlet-guide-vane position 
a. Thus, the fan model has no intrinsic dynamics. Inlet-guide-vane actuator 
dynamics are modeled as a second-order linear dynamical system. Fan-drive dynamics 
are not included in the present model because constant fan speed is assumed. 

The functional relations between fan-inlet and fan-outlet flows established by 
the Gumas fan model establish the relation between the outlet and inlet flows of the 
closed-return wind tunnel. This can be represented formally as a vector-function 
equation. In this solution, let V denote the vector of three independent flow 
properties 


V 


pw 

puoj 

e^po) 


(73) 


where subscripts 1 and N denote tunnel inlet and tunnel outlet, respectively. 
Also, let G be the vector function representing fan-outlet flow as a function of 
fan-inlet flow, fan speed, and inlet-guide-vane position. Then, the relation 
between (the vector of tunnel-inlet flow properties) and (the vector of 

tunnel-outlet flow properties) is expressed by 


Vi = G(V^,a,v) 


(74) 


— 1 

The relationships symbolized by vector function G and its inverse G determine 
the boundary conditions at the ends of the one-dimensional tube. Functions G and 
G””* are evaluated numerically by means of a computer subroutine. 

NUMERICAL SOLUTION OF FLOW EQUATIONS 

Equations (8), (16), and (18) are solved numerically to obtain a time-dependent 
solution of the wind-tunnel flow process. For notational convenience, these 
equations are written in conservation form as 


9t 


V(x, t) 



F(x,t) 


+ H(x,t) 


0 


(75) 


where 


puo) 


F = 


2 , 4 9u 

pu 0) + pco - - (iU 


B^puto + puo) - kto 


dT 

9x 


4 5u 

- 3 puco ^ 


(76) 
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and 


H 


-m' 



(77) 


A second-order explicit predictor-corrector method developed by MacCormack 
(ref. 14) is employed to obtain a numerical solution to equation (75). In this 
method, let the length of the tube be divided into N equal segments of length Ax 
and let At^ denote the time increment at the time step i. Let the diacritical 
marks ~ and denote predicted and corrected values, respectively, of the 
variables. Let subscript n denote the values of the variables at station x^^ 
where 


X = n Ax 
n 


(n — 0, . . . , N) 


(78) 


and let superscript k denote the values of variables at time t, where 




(k = 0, 1, 2, . . . ) 


i=1 


so that the predicted value is 


(79) 


^n = 


(80) 


and the corrected value is 


Vn = V(x^,t^) 


(81) 


MacCormack's method (ref. 14) proceeds as follows for n = 2, 3, ..., N-1; 


( 1) Predictor step; 


^+1 = - 


(2) Corrector step: 


yk+ 1 _ ^k 

n n 


(^n - ^-l)fc + 


At 




At 


(82) 


(83) 
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(3) Update step: 



( 84 ) 


Because partial derivatives occur in vector function F, additional differencing 
is necessary within and fJ^. In MacCorraack's formulation, these derivatives are 
written, respectively, as 


F^ = 

n 


p u to 
n n n 


[ p u^ + p - 4 P(u - u )-?— 1(0 
■^n n n 3 n+1 n AxJ n 

[ e^ pu +pu - k(T^.-T )-^ - 4 p(u - u )4— 1(0 
t,n n n n n n+1 n Ax 3 n+1 n AxJ n 


(85a) 


and 


~k 

F = 
n 


Pn^n^^ 


(85b) 


[pn’^n + Pn - I P<"n ’ ^n-1>7k]‘^n 
[«t,nPn^n + Pn% " ~ " I ’ ^-1>Er]"n 


The wind tunnel is a closed-return tube driven by a fan. As discussed previ- 
ously, the fan is represented by a vector-function relation G between vectors Vj^ 
and so that 


V, . G(V^) 


(86) 


From stations x^ to x^, the predictor-corrector steps are the following; 


(1) Predictor step; 


- 


' ^S-l)fe + At 


(87) 




r’ ■ 


( 88 ) 
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(2) Corrector step; 


^+1 = - 



(3) Update step; 



N 


G 


-1 



( 89 ) 


(90) 


(91) 


_ 1 

where G denotes the inverse of vector function G. 

The time increment At is recomputed after each time step with the relation 


0.9 Ax 

At = 

max ( I u I + a ) 


(92) 


where a is the local speed of sound and max denotes maximum value. 

Above a Mach number of 1.0, shocks form in the test section and produce rapid 
numerical changes which can lead to computational instability. Artificial viscosity 
is introduced into the conputation to smooth the shock discontinuity to a continuous 
change over several spatial increments. In the present treatment, an artificial- 
viscosity confutation method furnished by Turkel (ref. 15) is used. In this method, 
the following functions of time step k and spatial index n are defined as 


k _ S 

*^T Ax 


u,n 


k k 

u „ - u 
n+1 n 


V,n n+1 n 


= C^d’" 

V,n T u,n V,n 


The artificial-viscosity vector 



given by 


u^^-* = 


V,n 


(93) 


(94) 


(95) 


(96) 


(97) 


P'" 

V,n-1 
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is added to the updated state vector , given in equation (84), to yield the 

smoothed state vector 



n 


(n = 2, . . . , N-1) 


( 98 ) 


No smoothing is done at stations and Xjj- 

The aforementioned nimierical application of MacCormack's predictor-corrector 
method has been implemented on a CDC CYBER 203 digital coitputer. (See ref. 4 by 
Lambiotte and Howser. ) Vector equations (82) and (84) are well-suited for vectoriza- 
tion on this cortputer because each vector-addition and vector-multiplication opera- 
tion is performed by a single machine instruction for vectors of lengths up to 
64 000 elements. For the NTF model, the x-coordinate is divided into 512 segments. 

Because each of the three elements of vector for n = 2 through 511 in equa- 

tions (82) through (84) undergoes the same arithmetic operation, a single machine 
instruction is sufficient to perform the arithmetic operation on the entire set of 
1530 elements referenced in vectors through V 5 'ji» Separate operations are 

needed at stations 1 and 512 as given in equations (86) through (91). However, the 
elimination of 1530 separate addition or multiplication instructions for each vector 
addition or multiplication operation in equations (82) through (84) , along with loop- 
control instructions, results in a significant reduction in coirputing time. 

The proportion of scalar to vector operations in a program determines the extent 
of speedup over a serial program. Because the solution of the basic equations of 
flow is highly vectorized, a speedup factor of 50 to 60 is realized for the portion 
of the program corresponding to equations (82) through (84). However, because simu- 
lation of the actuator and feedback-control dynamics used in the present study is not 
vectorizable, inclusion of those features results in a larger proportion of nonvector 
operations. The conplete simulation, including actuators and controls, requires an 
average 0.006 sec of machine execution time per simulated time step At. The value 
of Ax is 0.295 m. From equation (92), it follows that 


At = 


0.266 

max ( I u I + a) 


(99) 


A typical value of At at a Mach number of 0.8 and total tenperature of 333 K is 
0.0004 sec. Thus, simulation of 1 sec of tunnel-control dynamics requires approxi- 
mately 15 sec on the central processing unit (CPU). At lower Mach numbers and tem- 
peratures, At becomes larger. For exanple, at a Mach number of 0.8 and tempera- 
ture of 167 K, At is 0.0005, and 12 CPU sec are required for each simulated tunnel 
second. In simulations of the Langley 0.3-Meter Transonic Cryogenic Tunnel, 256 seg- 
ments were used, which reduced execution times by a factor of 4 over that required 
for simulation of the NTF wind tunnel. 


INITIAL CONDITIONS 

Initial steady-flow conditions are obtained by integration of the steady-flow 
equations (A1) through (A3) of the appendix with an Adams method. (See ref. 16 by 
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Hindmarsh. ) Generally, some combination of steady-flow initial values of Mach num- 
ber, total tenperature, and total pressure at the test section are desired. However, 
values of state variables can be specified only at the tunnel inlet, along with 
GN 2 -vent flow rate Wq and LN 2 -in jection flow rate w^^. At steady flow it is seen 
that Wq must equal w^. Furthermore, the values of Wq and Wj^ must be such that 
fan work added to the enthalpies removed by LN 2 injection and GN 2 venting equals 
zero. Thus, values of tunnel-inlet state variables must be selected which will pro- 
duce a tunnel energy balance and the required values of M, T^, and p^ in the test 
section. 

It is convenient to employ as an indicator of energy balance at the fan. 

It can be shown that for steady flow, remains constant along the tunnel circuit 

except at the LN 2 spray bar and across the fan. Fan-outlet total tenperature is 
given by 



( 100 ) 


where 

rp fan conpression ratio 

T)^ fan efficiency 

fan-inlet total teitperature 

Because n is known and rp is conputed from the ratio of total pressures 


= P 4 . 1/P4. VI 

t y 1 t / N 


( 101 ) 


where p^ ^ denotes tunnel-outlet total pressure and p^ ^ denotes tunnel-inlet 
total pressure, p can be conputed from the integratecl solution of equations (A1) 

through (A3). Thus^ a procedure for obtaining initial conditions is as follows: 

(1) Select values of p^, u^ , e^ and w^; (2) integrate equations (A1) through 
(A3); and (3) conpute M,pg, p^ ,pg, ,j,g, and p. Since, at steady flow, p 

must equal 1 / by means of some type of organized search procedure, such as a 

gradient method, update p^, u-j, ®t '^G iterate until p equals 

and Hpg, Pt,TS' "^t,TS the required values. 

SIMULATION STUDIES 

The Langley 0.3-Meter Transonic Cryogenic Tunnel 

The Langley 0.3-Meter Transonic Cryogenic Tunnel (0.3-m TCT) has previously been 
simulated by means of a hybrid dynamic model by Thibodeaux and Balakrishna. (See 
ref. 17.) Although the scales and dynamic time constants of the NTF and 0.3-m TCT 
facilities differ by an order of magnitude, Mach number, tenperature, and pressure 
ranges are coitparable. Table III lists parameters for the 0.3-m TCT and the NTF. 
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TABLE III.- PARAMETERS FOR THE 0.3-m TCT AND THE NTF WIND TUNNEL 


Parameter 

Wind tunnel 

0.3-m TCT 

NTF 

. 2 

Test-section area, m 

0.124 

6.25 

Tunnel circuit length, m ... 

21.7 

151.2 

Tunnel volume, m^ 

14 

5597 

Plenum volume, m^ 

0.680 

907 

Mass of liner, kg 

3200 

465 000 

Temperature range, K 

80 to 325 

78 to 340 

Pressure range, atm 

1.1 to 6.0 

1.0 to 9.0 

Mach number range 

Maximum Reynolds number 

0.2 to 0.98 

0.2 to 1.2 

per 0.25 m 

Circuit time, sec; 

90 X 10^ 

120 X 10® 

M - 1, T - 300 K 

0.6 

4.8 

M - 0.2, T - 100 K 

3.0 
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Furthermore, the basic circuit geometries are similar, boundary interfaces are of the 
same type, and the fans, GN 2 vents, and LN 2 spray bars are in the same relative posi- 
tions in the 0.3-m TCT (fig. 7) and in the NTF wind tunnel (fig. 1). Because of 
reduced conputer execution times required for the smaller 0.3-m TCT, preliminary sim- 
ulation studies were performed using 0.3-m parameters and a simplified fan repre- 
sentation in the model. Some typical results of these studies are now described. 



Figure 7.- Schematic diagrcim of Langley 0.3-Meter Transonic Cryogenic Tunnel. 


Pulsed upsets of 3-sec duration in fan conpression ratio, GN 2 ~vent flow rate, 
and LN^-injection flow rate from steady-flow conditions have been simulated. Because 
fan drive and valve actuator dynamics were not included in this 0.3-m simulation, 
these upsets are represented as idealized rectangular pulses, which result in the 
flow responses shown in figures 8, 9, and 10. Computed responses to a 3-sec pulsed 
decrement in fan compression ratio appear in figure 8 at 277 K, 5 atm, and Mach 0.69. 
The computed responses, which show temporary decreases in total temperature, total 
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Time, sec 


Figure 10.- Computed response of 0.3-m TCT to LN„-flow upset 



pressure, and Mach number, are very nearly rectangular because of the omission of 
actuator dynamics. Figure 9 contains the computed responses to a 3-sec pulsed incre- 
ment in GN 2 -vent flow at 277 K, 5 atm, and Mach 0.7. 

Figure 10 shows the computed responses to a 3-sec pulsed increment in LN 2 flow 
at approximately 224 K, 2.1 atm, and Mach 0.67. These responses show an initial 
rapid decrease in total temperature followed by a slower recovery. Total pressures 
exhibit a very small initial drop followed by a sizable increase. Mach numbers show 
a pulsed increment during the LN 2 pulse followed by quick recovery to the original 
value. The computed temperature response exhibits temperature fronts which propagate 
at stream velocity as predicted theoretically. Temperature fronts were not observed 
e>q)erimentally in the 0.3-m TCT in tests which employed thermocouple temperature 
transducers with time constants of approximately 0.5 to 1.0 sec. Additional tests 
employing faster-response thermocouples are planned as of this writing. 

The computed temperature fronts decay after 3 to 4 tunnel circuit times because 
of heat transfer between the gas and the tunnel liner. Without gas-to-liner heat 
transfer, the amplitude of the steps would not be attenuated with time. 

Temperature fronts have been observed experimentally in the T'2 cryogenic wind 
tunnel at the Toulouse Research Center, Toulouse, France. (See ref. 2 by Blanchard, 
Dor, and Breil. ) The fronts seen at Toulouse, induced by a step change of -36 per- 
cent in LN 2 flow, were observed to propagate at stream velocity. They had decayed 
after two tunnel circuits, and the leading edges had become longitudinally diffuse 
near the end of their existence. Analytic solutions of the flow equations predict 
such temperature fronts, confirming the numerical result. Determination of their 
existence in the NTF must await experimentation after its completion. 


NTF Actuator, Sensor, and Control Studies 

The mechanical actuators for GN 2 ~vent valves, LN 2 -inlet valves, and inlet guide 
vanes in the NTF have been modeled by Gumas (ref. 13) as second-order linear systems. 
Valve stroke-flow calibrations have been provided as functions of local stagnation, 
supply, and vent pressures. Also, two types of high-accuracy pressure transducers 
were dynamically modeled by Gumas; a sonar mercury manometer and a quartz bourdon- 
tube transducer. Both systems have dynamic-response times of several seconds. Tem- 
perature transducers are modeled as first-order linear systems with time constants of 
0.63 sec. 

Three discrete sampled-data feedback control loops of the proportional-integral 
(PI) type with variable feedback gains for controlling stagnation pressure, stagna- 
tion tenperature, and Mach number in the NTF were developed by Gumas. Stagnation 
pressure is regulated by a GN 2 -vent -valve control loop; stagnation temperature, by 
an LN 2 -inlet-valve control loop; and Mach number, by a fan inlet-guide-vane control 
loop. Nominal sampling-time increments are 0.1 sec for the pressure and Mach number 
control loops and 0.04 sec for the temperature control loop. 

The Gumas simulations of actuators, transducers, and feedback controls have been 
appended to the distributed-parameter NTF model. In order for the computations to 
be synchronized, the ordinary differential equations describing actuator dynamics are 
solved with simple Euler integration for a time step At^^ equal to that of the 
MacCormack method (ref. 14) for solving the partial differential equations. The 
small value of At^ ensures adequate accuracy. 
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A large number of experimental simulation studies have been performed to eval- 
uate the performance of the Gumas control laws at selected values of Mach number, 
temperature, and pressure. Typical runs simulate programmed set-point changes under 
automatic control in Mach number, total temperature, and total pressure with the 
other two variables fixed. Performance criteria, which require minimum settling time 
for the controlled variable and minimum upset in the fixed variables, may thus be 
tested. Figures 11 through 14 show simulated results for selected test cases. Fig- 
ure 11 shows a controlled Mach number change from 0.8 to 0.9 at 2.04 atm and 
333.3 K. The Mach number set-point command occurs at 5 sec; the settling time to a 
Mach number of 0.900 is approximately 15 sec. Upsets in p^ at T^. are minimal, 
less than 0.0013 atm and 0.167 K, respectively. Conputed time histories of sensed 
test-section Mach number, total tenperature, total pressure, inlet-guide-vane 
position, liN 2 -valve position, and GN 2 -valve position are shown. 

Figures 12 and 13 illustrate a controlled temperature transition from 166.7 K to 
178.9 K for a Mach number of 0.8 at 2.04 atm. The temperature set-point command 
occurs at 5 sec, and the transition time is approximately 37 sec. Figure 12 contains 
actuator positions for guide-vane angle, LH 2 , and GN 2 , along with test-section Mach 
number, total pressure, and total temperature, the latter exhibiting well-defined 
temperature fronts. Pressure variation during the temperature transition is less 
than 0.0013 atm; Mach number variation is less than 0.004. Figure 13 shows the prop- 
agation of the temperature fronts, moving at stream velocity, at five stations along 
the NTF circuit. The fronts decay within six tunnel circuit times. 

Figure 14 shows a controlled pressure transition from 2.04 to 1.85 atm for a 
Mach number of 0.6 at 166.7 K. The pressure set-point command occurs at 5 sec. The 
pressure transition time is approximately 7 sec; however, Mach number and temperature 
disturbances require 11- and 13-sec settling times, respectively. It is seen that a 
significant temperature pulse of -3.3-K amplitude and 13-sec duration is introduced 
by the controlled pressure transition. Likewise, a Mach number upset of 0.010 ampli- 
tude occurs. Thus, although the Mach number and temperature control loops, which 
regulate inlet-guide-vane positions and LN 2 -inlet-valve position, do not appreciably 
interact with the GN 2 control loop, the latter control loop interacts significantly 
with the other two loops in regulating the vent-valve position. 

The simulated effects of a hypothetical guide-vane upset on steady controlled 
flow in the NTF are seen in figure 15. A 2.5° step pulse in inlet-guide-vane posi- 
tion of 2-sec duration is applied with Mach number controls fixed but temperature and 
pressure controls active. A significant Mach number transient results which decays 
after 10 sec. However, a train of temperature-impulsive disturbances of 2-K to 3-K 
amplitude is established which does not decay with time. It can be seen from the 
valve-position records that the reaction of the controls to the temperature impulses 
tends to reinforce and maintain their existence. This test case demonstrates the 
need for the carefully programmed Mach number transition control which was designed 
into the actual NTF control system. 

The effect of control-loop sampling rate on control performance was investigated 
in a series of simulations. The nominal values of the sampling intervals quoted 
previously were determined to be the largest values allowable for acceptable per- 
formance. Longer sampling intervals were found to cause slower control response, 
increased overshoot, and excessive oscillation. Sampling intervals of 0.5 sec or 
greater caused the control loops to become unstable. 

In a study by Armstrong and Tripp (ref. 18), multivariable-design techniques are 
applied to Mach number control of the NTF wind tunnel. In particular, optimal 


35 


<L> 

L. 

Z3 

J5 (o 

E 

a> 


cu 

o 


2.00 


J I I I L 


c 

o 


00 

1 

•4^ 

cn 

O) 


CD 

JQ 


E 


3 

C 


O 

aj 


S' 

“O 

oT 

cn 

c 

(O 

O) 

c 

(U 

> 

I 

o> 

;p 

*3 

O 


10 


-10 L 



I I I I I I 1 I I I I I I I 

0 4 8 12 16 20 24 28 


Time, sec 


Figure 11.- Simulation of controlled Mach number as predicted by NTF model 
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linear-regulator theory and eigenvalue-placement theory are employed to develop Mach 
number control laws, which are evaluated by using the distributed-parameter simula- 
tion described herein with actuator dynamics included. The resulting Mach number 
control law is found to reduce settling time significantly over that achieved by a 
conventional proportional-integral control loop. Details of the performance compari- 
sons described here appear in reference 18. 


CONCLUDING REMARKS 

The three-dimensional equations of fluid flow in integral form have been applied 
to a differential element of a one-dimensional tube of varying cross-sectional area. 
Heat-transfer effects due to thermal conductivity and viscous terms have been intro- 
duced to derive a one-dimensional system of partial differential equations analogous 
to the three-dimensional Navier-Stokes equations, which were applied to wind-tunnel 
geometry. The increase in entropy along the length of the tube was seen to be caused 
by a viscous momentum loss. This momentum loss, manifested as a loss in total pres- 
sure, was quantified from empirical pressure-loss relations obtained for diffusers, 
turning vanes, screens, and cooling coils in the National Transonic Facility (NTF) 
wind tunnel at the Langley Research Center. The cryogenic wind tunnels to be modeled 
contain a slotted test section whose analysis is based on a lumped model wherein the 
plenum is represented as a Iximped volume. The model is extended by separating the 
lumped plenum from the distributed test section and by distributing slot flow and the 
test-section length. 

The NTF and the Langley 0.3-Meter Transonic Cryogenic Tunnel are simulated by 
means of this model. A one-dimensional model of the fan adapted from a three- 
dimensional fan model is employed to relate the inlet-outlet values of mass, momen- 
tum, and energy flow rates in the one -dimensional tube. Pan compression ratio, 
determined by fan speed and an inlet-guide-vane system, establishes Mach number. 
Temperature control is effected by means of a liquid-nitrogen inlet spray bar and 
control valve. Pressure is maintained by a vent control valve. Inlet-guide-vane 
and control-valve actuators are modeled as second-order linear systems. 

The partial differential equations of the distributed parameter model are solved 
numerically by using an explicit, second-order, finite-difference, predictor- 
corrector method. The lumped ordinary differential equations describing actuator and 
control dynamics are solved by simple Euler integration synchronized with the finite- 
difference computation. 

The model has been employed in the development of multivariable control tech- 
niques based on optimal-regulator theory and eigenvalue-placement theory. ( See 
NASA TP-1887.) It has been employed extensively in the development of digital- 
process control algorithms for controlling Mach number, total pressure, and total 
temperature in the NTF wind tunnel. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
July 6, 1983 
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APPENDIX 


SOME STEADY-FLOW RELATIONSHIPS 

The differential equations for one-dimensional steady flow are 
^(puw) = m' 


(A1) 



pu (0 + pto) 


do) , . , 


(A2) 


d , , , 

— (e^puw + puw) = q' 

Equations (A1 ) and (A2) can be written, respectively, as 

du ^ u ^ _ _ u dw 
dx p dx 0 ) dx pw 

and 

du , 1 dp j ' - urn' 

y _ + _ ^ = 

dx p dx po) 

Total enthalpy h^ is defined as 


(A3) 


(A4) 


(A5) 



where e is the internal energy per unit mass. Equation (A3) becomes 
l^(h^puco) = q- 


or 

dht q' - h^m' 

u - — = 

dx po) 

For an ideal gas, the total enthalpy can be written as 


h 


t 




2 


u 


(A6) 


(A7) 


(AS) 


(A9) 
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With siibstitution of equation (A9), equation (A8) can be written as 


u d[£ _ u d£ Y ~ 1 pu du 

p dx p dx Y P 




- h^iti' 


pu 


(A10) 


This expression can be substituted into equation (A4) to yield 
du u_ ^ _ u ^ 2L ” ^ ^ ^ 2 


r(i' ■^1“'’”' -”i’) 


(All) 


dx YP YP<J^ 

Equations (A5) and (All) are solved simultaneously to eliminate du/dx so that 

(, . 1 . £!i! £}iL. * 1 (J,2) 

\ y p / dx udxu-^ YP<*>^ 2 ■'/ 

It is desired to obtain a steady-flow expression for the gradient of total pressure 
dp^/dx in terms of m* , j', and q' . Equation (A12) will be found to be useful to 
that end. The defining equation for p^ is 


/ 


(A13) 


where 




(A14) 


Differentiation of equations (A2) and (A13) and considerable manipulation yield 


dp 


t 


dx 



I pp ^ 

Y p dx 


(A15) 


One may substitute equation (A12) into equation (A15) 
that 

1 

^ = + j, JBh!) (2 + JLJlJ. _ I2Ll . JL^ 

dx 2 \ p / |_\u Y P“ / w Y 



1 Pu^V p 

Y P 


so 


(A16) 
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Finally , note that 


= Y - .I „2 


(A17) 


so that equation (A16) becomes 


dx 




(A18) 


Thus, dp^/dx is a linear function of m' , j', and q' . 

Now, obtain an expression for the gradient of entropy ds/dx for steady flow. 
Differentiate equation (A6) to obtain 


dh , , , , 

t de . dv , dp , du 

- — = ':^+Pi;~+v^+u — 
dx dx dx dx dx 


(A19) 


By the second law of thermodynamics, the differentials de, ds, and dv for a pure 
substance are related as 


de = T ds - p dv 


(A20) 


Substitution of equation (A20) into equation (A19) gives 


t „ ds dp , du 

— + T — + v^+u 
dx dx dx dx 


(A21) 


Expressions for dh^/dx and u ^ ^ obtained from equations (A5) and (A9). 

Substitution of these expressions into equation (A21) gives 


ds 

dx 


1 | um' - i ' 

T y pu 



Equation (A22) may be rewritten as 


ds 

dx 



- uj- 


+ q* 


(A22) 


(A23) 
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